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Abstract. We consider the Ginzburg-Landau functional defined over a bounded and smooth 
three dimensional domain. Supposing that the magnetic field is comparable with the second 
critical field and that the Ginzburg-Landau parameter is large, we determine a sharp asymptotic 
estimate of the minimizing energy. In particular, this shows how bulk superconductivity de- 
' creases in average as the applied magnetic field approaches the second critical field from below. 

_ , Other estimates are also obtained which allow us to obtain, in a subsequent paper [f 8| . a fine 

^ ' characterization of the second critical field. The approach relies on a careful analysis of several 

, limiting energies, which is of independent interest. 

\o 

<, 

r~| ! 1.1. Introduction. Since the publication of the pioneering books [3 [T2], the analysis of the 

\ Ginzburg-Landau functional of superconductivity became the subject of a vast mathematical 

literature. In particular, in the presence of an applied magnetic field and in an asymptotic 
limit where the Ginzburg-Landau parameter tends to oo, different advanced mathematical tools 
indicates three critical values of the applied magnetic field, usually denoted by Hqx , and Hc-j ■ 
[ These critical fields may be described in a rough way as follows (see e.g. the book of de Gennes 

' |10)). If a superconducting sample is subject to a constant applied magnetic field of intensity 
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CN ' H, then as long as H < Hc^^ the sample is in a pure superconducting state and repels the 

magnetic field. If the field is slightly increased above Hc^ , the sample looses superconductivity 
I in point defects called vortices, whose number increases as long as the field H is increased all 

the way up to the critical value Hci', this phase of the superconducting sample is the mixed 
phase. Increasing the field above destroys superconductivity in the bulk of the sample, but 
as long as the field H is below Hc-^, the sample carries superconductivity along the surface. 
^ I Finally, if the field is increased past Hq^, superconductivity is lost in the sample which switches 

^ I to the normal state. However, to establish these results in a rigorous mathematical framework, a 

variety of notations and methods are required. We mention here that the analysis of the presence 
of vortices is the subject of the monograph |25| and references therein, while that of surface 
superconductivity is that of |12) (and references therein). The methods originally developed to 
the analysis of the Ginzburg-Landau functional are also appearing relevant in the analysis of other 
mathematical models of condensed matter physics like superfluidity, Bose-Einstein condensates 
and liquid crystals (see [Il[9l[23] and the references therein). 

In this paper, we focus on the analysis of the Ginzburg-Landau functional in a three dimen- 
sional domain and as the Ginzburg-Landau parameter tends to oo. The regime of the applied 
magnetic field considered here is when the field intensity varies close to and below the second 
critical field ■ In two other subsequent papers |18| [20] , the analysis is completed to cover 
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variations of the magnetic field intensity in the regime ^ Hc^ and increasing all the way up 
to Hc^- The novel issue in the results is that previous ones where only present for the two 
dimensional Ginzburg-Landau functional, see \12\ I25|. 

Many questions which are answered successfully for the two-dimensional Ginzburg-Landau 
functional remain open for the three-dimensional functional. This includes, among other things, 
a determination of the first critical field Hc^ as precise as the one present in two-dimensionsQ 
Among the other open questions, we mention the analysis of the mixed phase and the analysis 
of the critical field (that is part of the subject of this and the subsequent papers |18| I20|). 

However, in three-dimensions, the analysis of the third critical field Hc^, already started in 
|21| , then a sharp characterization of it is given in |14| . The regime of an external magnetic field 
strength close to and below Hco, is discussed in the papers [5]|22]. In particular, it is proved that 
superconductivity is confined to the surface. These results will be significantly improved in the 
subsequent paper [TS]. 

In [5] , interesting estimates are proved in the regime of applied fields close to and below , 
which also is the regime we treat in this paper. We improve the estimates in [5] by proving sharp 
estimates. The approach we follow is considerably different from [5] and mainly variational in 
nature. 

1.2. The functional and main results. We consider a bounded and open set C M'^ with 
smooth boundary. We suppose that Q models a superconducting sample subject to an applied 
external magnetic field. The energy of the sample is given by the Ginzburg-Landau functional. 



Jn 



|(V-m/7A)V'p-KVl' + yl 



dx 



+ k^H'^ I |curlA-/3|2(i2;. (1.1) 

Jr3 



Here k and H are two positive parameters; k (the Ginzburg-Landau constant) is a material 
parameter and H measures the intensity of the applied magnetic field. The wave function (order 
parameter) € H^{Q;C) describes the superconducting properties of the material. The induced 
magnetic field is KiifcurlA, where the potential A G i7jj^p(R^), with i7jjyp(R^) being the 
natural variational space for the problem defined in (|1.2p below. Finally, f3 is the profile and 
direction of the external magnetic field that we choose constant, /3 = (0,0, 1). 

The space ij^(IR^) is the homogeneous Sobolev space, i.e. the closure of C^(M^) under the 
norm n i— )• ||^t||/j-i(ig3) := ||Vit||/^2(]R3). Let further F{x) = (— X2/2, xi/2, 0). Clearly divF = 0. 

We define the space, 

^div,F(lR^) = {A : divA = 0, and A-F£H^{R^)}. (1.2) 

Critical points {ip,A) G H^{n;C) x ij^j^p(]R'^) of are weak solutions of the Ginzburg- 
Landau equations, 

-{V -iKHAfip = k'^{1-\iP\'^)tP in 

curl2A = —ImN^ (V -iKH AU)ln in (13) 

kH 

u-(V - iKHA)ij = on dfl, 

where 1q is the characteristic function of the domain Q, and 1/ is the unit interior normal vector 

of dn. 



^It is only recently that results related to this regime of the magnetic field and valid for general domains 
appeared in [6]. Earlier contributions include a candidate for the expression of Hci in the case of the ball [3], and 
an expression of Hci in the case of 'thin' shell domains [S]. 
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Our main results include a sharp asymptotic estimate of the ground state energy which is 
defined as follows, 

Eg.st(K,^) =inf{f3°(V',A) : {i;, A) € H\n;C) x HI,.,f(.^^)} ■ (1.4) 

The asymptotic estimate of Eg^st{K, H) involves an auxiliary function g : [0, oo) — >■ [—^,0] intro- 
duced in (|2.5p below. The function g is increasing, continuous and g{b) = for all 6 > 1. 

Theorem 1.1. Let < Amin < Amax- There exist positive constants C and kq such that if 

K 

then the ground state energy in (|1.4p satisfies, 

Eg.st(K, H) - g{H/K)\n\K^ < Ck^I'^ . (1.5) 

Theorem 11.11 is an extension to three-dimensional domains of the analogous Theorem 1.4 in 
established for two-dimensional domains. Let us emphasize that the limiting function gi^-) 
in Theorem 11.11 is the same as that in [26], up to an additive constant. 

In the particular regime when U = «;-|-o(k). Theorem 11.11 tells us that Eg.st('t, -f^) = o(k^) as 
K — 7- oo. In this case, the leading order term in the ground state energy is given in the subsequent 
part of this paper |18j . Actually, we will prove that there exist two constants -Esurf < and £'2 < 
such that, 

Eg,st(K,i?) =£^surf'« + £^2|^^|[K-i^]+ + o(^max(K,[K-F]^)y (1.6) 

Here we use the function = max(x,0), x G M. We stress that the constant i?surf depends 
(only) on the boundary of 0, while the constant Ei is universal. Part of the necessary estimates 
to prove ()1.6p are given in the last section of this paper. 

Furthermore, in the subsequent part of this paper |18) . we determine leading order estimates 
of the ground state energy when the magnetic field increases all the way up to the critical 
value i^c'g . This gives an overall understanding of the ground state energy of three dimensional 
superconductors subject to strong magnetic fields, with a precise description of the transition 
from hulk to surface regime as well. 

The next theorem concerns the behavior of order parameters for critical configurations. We 
obtain a sharp asymptotic estimate on the average of the order parameter on small cubes. 

Theorem 1.2. Suppose the assumptions of Theorem \l.l\ are satisfied. There exist positive con- 
stants C , Rq and kq such that, if 

Kf) < K , RqKT^ < i < - , 

and Qi is a cube of side-length I and parallel to the external magnetic field (3, satisfying Qi C Vt, 
then the following is true. 

(1) If{ip,A) G i?i(f];C) X Hl-^^piM.^) is a solution of ([OD, then, 

^ ' ' "^dx<-2g{H/K) + c(K-^/^ + ^] , as k^oo. (1.7) 



\Qi\ jQi V 

(2) If{^p,A) G H^{n;C) X if]j^p(M3) is a minimizer of ([LU, then, 



\Qi 



\'4j\'^dx + 2g{H/K) 



Qt 



<C{k'^^^ + ^], as (1.8) 
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Remark 1.3. Theorem 11.21 still holds true (with a different error term) if cubes are replaced by 
more general domains that can be approximated by cubes of side- length i and <C ^ <C 1 as 
K — 7- oo. In particular, it is true for balls of radius i. 

Remark 1.4. The function g is a continuous increasing function, and g(b) = for all 6 > 1 (see 
Section [2. ip . Furthermore, there exists a universal constant a G (0, ^) such that, 

V6G [0,1], a(6-l)2 < \g{b)\ < lib - if ■ 

Theorem 11.21 then shows that as the applied magnetic field approaches (i-e. H ^ bn and 
6 — )• 1_), superconductivity decreases in the bulk like \/l ~ b. 

Remark 1.5. In connection with the result of Theorem 11.21 Almog [5] proves that if the magnetic 
field is such that H = bn + o{k) with 6 < 1, and if (ip, A) is a solution of the Ginzburg-Landau 
equations ()1.3p . then 



/ \^l^\Ux<- + 2 
Jn « 



\ft\ 



(1,9) 



The estimate ()1.7p we give in Theorem 11.21 is an improvement of (jl.Op in the sense that it 
estimates the average of \ip\^ on small cubes and to leading order. 

In the case b = 1, the estimate of the ground state energy in ()1.6p yields an estimate stronger 
than the ones in (jl.9p and Theorem 11.21 We refer to |18| for the precise statement. 

Remark 1.6. In light of the result of Theorem 11.21 it is natural to expect that for any solution 
(V', A) of the Ginzburg-Landau equation ()1.3p . the order parameter ip satisfies 

1/2 



l^{lu^) < Cmax 



where C is a constant, S € (0, 1) and oj,^ = {x £ ft : dist(x,5r2) > k~^}. Bounds of this type 
were obtained for the two dimensional equations in |13| [T6| 117) . but we are not at present able 
to obtain them in 3D. 



The general technique used in this paper is mainly based on the methods of |12) . After deter- 
mining a priori estimates of the solutions, we reduce the problem, via localization techniques, to 
that of a three dimensional model problem. The latter is linked to a two-dimensional problem, 
which has been studied in [21 [26] . 

The main concern in [U [26] was the leading order term in the ground state energy for the two- 
dimensional model problem, without attention to estimates of the remainder terms. However, 
precise estimates of those remainder terms are necessary in order to obtain the correct leading 
order term in the energy (|1.4p . 

We perform a very careful analysis of the two dimensional limiting problem thereby obtaining 
estimates that, on the one hand, are stronger than those appearing in [H [26], and on the other 
hand, when inserted in the original functional in (jl.ip . their contribution to the remainder terms 
is negligeable compared with the expected leading order term in the energy. On a technical 
level, a key tool in the control of remainder terms is an L°°-bound for solutions of a reduced 
Ginzburg-Landau equation in the plane, which was obtained in |13) . 

Let us point out that, using less sophisticated tools than the ones used in [2] [26], we recover 
in a unified approach all the results obtained in [2] concerning the limiting problem. It seems 
that our approach works for limiting energies arising in other contexts. In |18| . through the same 
approach, we identify the limiting problem of the surface regime. 

The paper is organized as follows. In Section [21 we analyze the limiting energy. In Section [3l 
we give asymptotic a priori estimates on solutions of the Ginzburg-Landau equation (jl.3p . In 



3D GINZBURG-LANDAU FUNCTIONAL 



5 



Section [U we prove local energy estimates in small cubes for minimizers of the energy in (ll.ip . 
Section [5] concludes with the proofs of the theorems announced in the introduction. Section [U] is 
devoted to additional estimates related to the case when H is close to k, which will be used in 

m- 

Remark on notation: 

• Throughout the paper, we write £ for the functional in (jl.ip . 

• The letter C denotes a positive constant that is independent of the parameters k, and H, 
and whose value may change from line to line. 

• If a{K) and 6(k) are two positive functions, we write a{K) <^ if a{K)/b{K.) — )■ as 

K —7- OO. 

• If a(K) and 6(k) are two functions with ^ 0, we write a{K) ~ if a{K.)/b{K) — )• 1 
as K — OO. 

• If a(n) and are two positive functions, we write a{K) ~ if there exist positive 
constants ci, C2 and kq such that cib{n) < a(n) < C26(k) for all k > kq. 

• If X G M, we let = max(a;,0). 

2. The limiting energy 

This section contains the study of the large area/volune limit of Ginzburg-Landau functionals 
with constant magnetic field in 2 and 3 dimensions. These models are fundamental for the results 
of the paper. In the 2D case we compare different boundary conditions (Dirichlet in Section 12.11 
and periodic in Section I2.1.2P and introduce the Abrikosov energy corresponding to restricting 
the functional to periodic functions in the lowest Landau band in Section 12.1.31 Finally, in 
Section [2.21 we reduce the 3 dimensional case to the 2D one. 

2.1. Two-dimensional limiting energy. 

2.1.1. Reduced Ginzburg-Landau functional and thermodynamic limit. 

Given a constant b > and an open set P C M^, we define the following Ginzburg-Landau 
energy, 

Gb.v{u) = j (^&|(V-iAo)n|2-|n|2 + l|n|^^ dx. (2.1) 

Here Aq is the canonical magnetic potential, 

Aq{xi,X2) = ]^{-X2,xi) , V X = (xi,a;2) G . (2.2) 

We will consider the functional Gi,^x> first with Dirichlet and later with (magnetic) periodic 
boundary conditions. It will be clear from the context what is meant. 

Consider the functional with Dirichlet boundary conditions and for 6 > 0. If the domain T> 
is bounded, completing the square in the expression of G^p shows that Gh^v is bounded from 
below. Thus, starting from a minimizing sequence, it is easy to check that Gh^v has a minimizer. 
A standard application of the maximum principle shows that, if u is any minimizer of G;,^-pj then 

\u\ < 1, in V, (2.3) 

see e.g. [25] . 

Given i? > 0, we denote by Kr = {-R/2, R/2) x {-R/2, R/2) a square of side length R. Let, 

mo(6,i?)= inf Gb,KM- (2.4) 

u&HI(Kr;C) 

The main concern of this section is the proof of the following theorem. 
Theorem 2.1. Let mo{b,R) be as defined in ()2.4p . 
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(1) For all b > 1 and R > 0, we have mo{b, R) = 0. 

(2) Ifb€z [0,oo), there exists a constant g{b) < such that, 

(3) The function [0, +oo) 9 6 i— )• g(b) is continuous, non- decreasing, concave and its range is 
the interval [—^,0]. 

(4) There exists a universal constant a G (0, ^) such that, 

V6g[0,1], aib-lf <\g{b)\<^{b-l)\ (2.6) 

(5) There exist universal constants C and Rq such that, 

y R>Ro,ybe[0,l], g(^) < "^o(fe,i^) <ff(fe) + |. (2.7) 

The major part of Theorem 12.11 is obtained by Sandier- Serf aty [26] and Aftahon-Serfaty [21 
Lemma 2.4]. However, the estimate in (j2.7p is new. We give an alternative proof of Theorem 12. II 
through a direct approach^. 

An important key-ingredient in the proof of Theorem 12.11 is the following abstract lemma. 

Lemma 2.2. Consider a decreasing function d : (0, oo) — )■ (— oo,0] such that the function 
f : {0,oo) 3 £^ ^ eM. is bounded. 

(1) Suppose that there exist constants C > and > such that the estimate 

fini) > /((I + a)i) -C(^a + -1^ ) , (2.8) 

holds true for all a G (0, 1), n G N and i > io- 

Then f{l) has a limit A as t ^ oo. Furthermore, for all I > IIq, the following estimate 
holds true, 

Of 

m<A + ^,. (2.9) 

(2) Suppose that there exist constants C > and > such that the estimate 

f{n£)>fi£)-j, (2.10) 

holds true for all n and £ > Iq. 

Then f{£) has a limit A as £ ^ oo. Furthermore, for all £ > 2£q, the following estimate 
holds true, 

on 

fi£)<A + — . (2.11) 

Proof. We prove the first conclusion of the lemma. As the proof of the second conclusion is 
similar, we do not write the details. 

We denote by j4 = limsup/(^). We know that A is finite since the function / is bounded. 

£~>-oo 

Also, A < since d is non-positive. We will prove that liminf f{£) > A. This will give us that 

f{£) has limit A as £ ^ oo. 

Suppose that e G (0, 1) is a given constant. We may select £o > l/(e^) such that f{£o) > A-e, 
and the estimate in ()2.8p holds true for all £ > £o. 



^In [2] [26], the authors link the ground state energy mo (6, R) to that of a full Ginzburg-Landau energy (where 
the magnetic field is an unknown), then for the latter, they analyze a simultaneous limit of both large area and 
large Ginzburg-Landau parameter. 
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For each £ G (£o/(l+£), oo), let n G N be the unique natural number satisfying n£o < {l+e)i < 
(n + l)£o- By monotonicity of ^ i-)- d{£), we have, d{£) > d{{n + 1)£q/{1 + e)). By definition of 
f{£) we get, 



(n + 1) 



'0 



/ 



(1 + 4 

Using the bound in (I2.8p with a = e, we deduce that. 



[n+l) 



'0 



1 + e 



(n + 1)4 



(1 + eK 

We notice that by our assumption on n, 



fi£o) -C{e + 



{n + l)£o 
{l + e)£ 



< 



n£(] 



< 1 + 



{l + e)£ 
2 + £o 



+ 



(1 + ^) 



2nio 

{l + e)£ ' {l + e)£ 



+ 



Since 



< 0, we deduce that, 

m > ( 1 + 



2 + . 



fi£o) -C(e + 



eHl 



Consequently, we get by taking liminf on both sides above, 



Taking e — )■ 0+, we get that 



liminf /(£) >A-C{e + e^ 



liminf /(^) >A. 



The estimate in ()2.8p applied with arbitrary n G N and a G (0, 1) yields, 

1 



f{nL)>f{{l + a)L)-C{a + 



for all L > £q. Making n — t- oo we get. 



x2L2 



A> f{{l + a)L)-C\^a + 
Selecting L = £/{\ + a) and a = £~'^l'^ we obtain that. 

This finishes the proof of assertion (1) in Lemma [2.21 

In the next lemma, we give rough bounds on the energy mQ{b, R). 

Lemma 2.3. For all b > and R > 0, we have, 

11 - b]lR' 



□ 



< mo{b,R) < 0, 

where mQ{b,R) is defined in (|2.4p . 

Furthermore, there exist universal constants a G (0, ^) and C > such that, for all b G [0, 1) 
and n G N, we have, 

(2.12) 



mo(6,2V2^n) ^2 ^ " 

■ < — a(l — o) + 



87rn2 



n 
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Proof. Let u E Hq{Kji). We extend u to a function u G H^{R'^) by setting u = in \ Kr. 
Then, 



dx . 



Recall that the lowest lowest eigenvalue of the Schrodinger operator with constant magnetic field 
— (V — iAo)^ in L^(M^) is equal to 1. So, we get by the variational principle, 

G,,K«M>^^((6-1)|SP + ^|S|^) dx = (^(6-l)|t.|2 + i|n|4) dx . 

Since u G Hq{Kji) was arbitrary this clearly gives the lower bound of the lemma when 5 > 1. 
When < 6 < 1, we can complete the square to get, 

((6-l)|np + l|u|4) rfx = i^_ ^^(b-l)-\u\^f-{l-bf)dx 

>-\{l-hfR\ 

which thereby finishes the proof of the lower bound of the lemma. 

Since Gh^Kni^) = 0, we get the trivial upper bound mo{b, R) < 0. This proves the first 
estimate in the lemma. 

The proof of the estimate in (j2.12p consists of computing the energy of a test function con- 
structed in |26| . A similar bound (more in line with the techniques of the present paper) can be 
obtained by using the 'lowest Landau band' functions introduced below. We will briefly sketch 
such a calculation. With R = 2\p2/K the space L/j defined in Proposition 12.91 is non-empty and 
contains a non-zero magnetic periodic function v € -£'/j=2v^ ' '^^^'^^ the space Er is defined 



in (|2.20p . By defining u = X\/l — bv/\\v\\i^2(^XR)j '^^ may select A > sufficiently small and a 
constant a G (0, 1/2) such that, 

Gb,K^iu) = -a{l - bf < 0. (2.13) 

We can take the magnetic periodic function u and cut it down to a box of size nR (in order to 
satisfy the Dirchlet boundary) condition. Upon calculating the energy on KnR of this function 
we get the bound (|2.12p . where the last term comes from the localization error. 

This finishes the proof of the lemma. □ 

The next lemma establishes monotonicity properties of mQ{b,R) with respect to R. 

Lemma 2.4. For all b > 0, the function R i— t- mQ(b,R) defined in (j2.4p is monotone non- 
increasing. 

Furthermore, there exist universal constants C > and £q > such that, for all b > and 
i > io, we have the estimate, 

mQ{b,nR) ^ mo{b,R) C{b + 1) 



Proof. We start by proving that the function R i— )• mQ{b,R) is monotone. Let r > and 
u S Hq{Kr). We extend u to a function u G //q (i^'/j+j.) by setting n = outside Kr. In this 
way, we get, 

Gh,Ka{u) = Gb,KR+ri^) ^ '^o{b, R + r) . 
Since u G Hq{Kr) was arbitrary this proves the montonicity of mQ{b,R) with respect to R. 
We prove the lower bound of the lemma. If j = (ji, 72) G we denote by 
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where 



' 2m + 1 — n 1 2m + 1 — n 1 
VmGZ, /„=( + - 



For all i? > 0, we set 

Let J = {j = (ji, js) G : < ii,j2 < n - 1} and KnR = {-nR,nR) x {-nR,nR). Then 
the family {ICRj)j^j is a covering of KnR, and is formed exactly of n? squares. Let UnR be a 
minimizer of Gb^K^R, i-^- Gii^k^rWur) = mQ{h,nR). We have the obvious decomposition, 

/ \UnR\^dx = ^ \UnR\'^dx. (2.14) 



Let xrj- E C^(M^) be a function satisfying, 

0<Xi?j<l, inR^, supp XRj- C /Ci?j , XR,i = 1 in -fCR_i,j , |VxR,j| < C in , 

for some universal constant C. 

Since satisfies the equation — 6(V — iAo)^Unfi = (1 ~ \unR\'^)unR, it results from an 

integration by parts that Gb^Knui'^nR) = / WnR\^ dx and. 



\0\\y - lJ\Q}XR,jUnR\- - \XR,jUnR\- +XR,j\UnR\') CX = j \\\/ XR,j\UnR^ 

Using that \unR\ < 1, it is easy to check that (with a new constant C) 



/ - i-^o)XR,jUnR\'^ - \XR,jUnR\^ + XR,j\UnR\'^) dx = b ||VxR,j|n„/ 



dx . 



Gb,K„R{XR,jUnR) < -\ [ \UnR\^dx + G{b+l)R. 



(2.15) 



Since each XRJ'^nR has support in a square of side length R, we get by magnetic translation 



invariance, 



Gb,K„RiXR,jUnR) > mo{b,R) . 
We insert this lower bound into ()2.15p . then we take the sum over all j G J. In this way we get. 



niQ 



= -\ I \unR\^dx + Cn'{b+l)R. 
To finish the proof of Lemma 12.41 we just recall that, 

mQ{b,nR) = Gb,K^R{unR) = -]: [ \unR\'^ dx . 



The next lemma treats the specific case 6 = 0. 

^2 

Lemma 2.5. If b = and R> 0, then ■mQ{b,R) = — —. 

Proof. Let u G HI{Kr). It is easy to see, by completing the square, that if 6 = 0, 



□ 



G,,Kr{u)= I l{\u\'-l)'dx-^. (2.16) 
JKr 2 2 
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Consequently, we get mo{b,R) > when & = 0. We will prove that mQ{b,R) < —R^/2. If 

n > i? is a natural number, let u„ E C^^Kpt) be a function satisfying, 

< u„ < 1, ii„ = 1 on Kj^ i. 

Then moib, R) < Gh^Kniun) for all n > R. It is easy to check that Lebesgue's dominated conver- 
gence theorem gives, lim GhKaiun) = > thereby proving that mo(6 = 0,R) < —R^/2. □ 

Proof of Theorem ] 2. 1\ The proof is separated into several steps. 

Step 1. Notice that the first conclusion in the theorem is a direct consequence of Lemma [2.31 
7Tl()(b R^ 

We prove that ^ has a limit g(b) as i? — )■ oo. Thanks to Lemma [2.51 we have nothing to 

prove when 6=0. Also, Lemma [2.51 gives us that ^(0) = —1/2. 
Suppose that < 6 < 1. Let 

/6(i?) = ^^^4^' MR) = ^o{b,R). 

Thanks to Lemmas 12.3! and I2.4| the functions and di, satisfy the assumptions of Lemma 12.21 
Therefore, we conclude that fb{R) has a limit g{b) as R ^ oo. Furthermore, since b £ [0, 1], we 
infer from the estimate of Lemma I2.2| 

f,{R) < g{b) + ^ , 

where C is a universal constant. This proves the upper bound in assertion (5) of Theorem 12.11 
Step 2. In this step, we prove the lower bound rriQ^b, R) > g{b)R? . Let n G N and 72 > 0. 
We use the notation Kr = (-f , f ) x (-f , f ) and Kr = {0,R) x {0,R). Let u £ HI{Kr) be 
a minimizer of G^j^^. By magnetic translation invariance, we get mQ(b,R) = Gj^j^^{u). We 

extend u to a function u G HQ^KnR) by 'magnetic periodicity' as follows, 

u{xi + R, X2) = e'^'-'^/'^u{xi,X2) , u{xi,X2 + R) = e~^^''^l'^u{xx,X2) ■ 
Using magnetic translation invariance, it is easy to check that 

^b,K„a(^) = n^mo{b,R) . 

Consequently, we get mo{b,nR) < n?mQ{b, R). We divide both sides of this inequality by n^i?^ 
then we take the limit as n — )• 00. That gives us g{b) < mQ{b, R). 
Step 3. 

In this step we prove that the function g is increasing. Let 6 > and e > 0. Since e > 0, it 
is easy to see that mo{b + e,R) > mo{b,R). Dividing both sides of this inequality by R^ then 
taking i? — )• 00 we get g{b + e) > g{b). 

Step 4- 

In this step we prove that g is concave and continuous. The concavity of g is straight forward. 
Upon writing 

moib,R)=mf j (^6|(V - iAo)up - |m|2 + i|u|^^ dx, (2.17) 

we see that b 1— t- mQ(b,R) is the infimum of a family of affine functions and therefore concave. 
So g{b) is the pointwise limit of concave functions and therefore concave. 

The concavity of g implies continuity except at the endpoint 6=0. Since ^(0) = — 1/2 and 
g is non-decreasing it suffices to prove that lim sup (7(e) < —1/2. Consider e > 0. In Step 2, we 
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proved that 

for all R>1. Let tt E H^{Kr). From (f2J8]) we get that 
lim sup (7(e) < i?~^limsup / | e|(V — lA 



mo(g,i?) 

9(e) < — ^2 — ' (2.18) 



= R~^J (^-H^ + dx. (2.19) 

Since n is arbitrary the result follows from Lemma [2.51 

Step 5. The assertion (4) in Theorem 12.11 follows directly from Lemma [2.31 

This finishes the proof of Theorem 12.11 □ 

2.1.2. Periodic minimizers. 

RecaU that for each R > 0, Kr = {-R/2,R/2) x {-R/2,R/2) is a square of side length R. 
We introduce the following space, 

^i? = |u G i7/„,(IR2;C) : u{xi + R,X2) = e'^^'/Mxi,X2) 

u{xi,X2 + R) = e-'^'''/^u{xi,X2) , V (xi, xa) G . (2.20) 



Notice that the periodicity conditions in (|2.20p are constructed in such a manner that, for any 
function u £ Er, the functions \u\, |Vao^I and the vector field uVaq^ are periodic with respect 
to the given lattice. 

Recall the functional v in (|2.ip above. We introduce the ground state energy, 



mp(6,i?) = inf Gi,,kM. (2.21) 

Starting from a minimizing sequence, it is easy to see that when 6 > 0, G^^Kj^ admit minimizers 
in the space Er. Various properties of the minimizers are collected in next proposition. 

Proposition 2.6. Let b > and let u G Ejj be a minimizer of the functional Gb^Kn, ^-e. 
Gb,Kii{u) = mp{b, R). Then u has the following properties. 

(1) u is a solution of the equation, 

- 6(V - iAo)\ = (1 - |n|2)u, inM^ (2.22) 

(2) There exists a universal constant Cmax > such that 



u||L->m2) < min l,Cn 



1 



1/2N 



(2.23) 



Proof. Since u is a minimizer of Gb^Kj^, then u satisfies the equation in (|2.22p in Kji. Using the 
periodicity of u together with the explicit definition of the magnetic potential Aq, it is easy to 
check that u satisfies the equation in M^. 

Since |n| is periodic, then ||ii||^oo(]g2) = |n(xo)| for some xq G Kjj. It results from a standard 
application of the strong maximum principle that |u(a;o)| < 1- This proves that ||ii||2,oo(iR2) < 1. 

Since u is a bounded solution of (|2.22p . then Property (2) of the proposition is a straight 
forward application of Theorem 3.1 in |13| . □ 

In the next proposition, we exhibit the relation between the ground state energies mo(6, i?) 
and mp{b, R). 
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Proposition 2.7. Let mQ{b,R) and mp{b,R) be as introduced in ()2.4p and ()2.2ip respectively. 
For all b > and R > 0, we have, 

mQ{b, R) > mp{b, R) . 

Furthermore, there exist universal constants eo G (0, 1) and C > such that, if b > 1 — eo and 
R>2, then, 

mo{b,R) <mp{b,R) + C[l-b]+R. (2.24) 

Proof. Since every function u £ Hq{Kji) can be extended by magnetic periodicity to a function 
in the space Eji, we get immediately that mo(6, i?) > mp{b,R). 

We prove the upper bound. Suppose that R> 2. Let XR £ C^{M?) be a function satisfying, 

0<Xij<l inM2, snppxRC Kr, xr = in AVi , I^Xr} < M mR\ 

for some universal constant M. 

Let u G Er be a minimizer of Gh^Kn^ i-e- Gh^Kniu) = mp(b,R). Let / = XR'^- Then 
/ G Hq{Kji) and consequently 

mo(6,i?) <G;,,,^«(/). 
The equation for u and an integration by parts yield, 



G,,KAf)= I b\\VxR\n\Ux + \ [ \xRu\Ux- [ Xr\ 
mp(b,R) = Gb,Kn{u) = -\ I 



u\ dx . 



and ^ ^ 

\u\^ dx . 

Using the properties of XR — and the bound ||ti||oo ^ 1 — we deduce the following upper bound, 

mo(6, R) < mp{b, R) + A{M^ + l)R\\u\\l^^Ka) ■ 
We get (|2.24p by inserting (|2.23p in this estimate. □ 

The next proposition gives a uniform upper bound of mp{b, R) which is interesting when b is 
close to 1 and R is large. It is a key-ingredient in the proof of Theorem 12.121 below. Actually, 
this bound is also true with mp{b, R) replaced by mo(6, R) but we will not need this version. 

Proposition 2.8. Let g[b) andmp{b,R) be as defined in (j2.5p and (|2.2ip respectively. There 
exist universal constants cq £ (0,1), Rq > 2 and C > such that, if 1 — e^ < b < 1 and R> Rq, 
then, 

mpib,R) ^ 5)3/2 
R^ ^ 

Proof. Let n > 2 be a natural number. If a G (0, 1) and j = (ji, j2) G Z^, we denote by 

Ka,j = Iji X , 

where 

«, „ ^ _/2m+l-n (1 + a) 2m + 1 - n (1 + a) 
VmGZ, Im-[^ - —, - 

Consider a partition of unity (xj) of such that: 

X]lXil^ = l, 0<Xi<l inM^, suppxjC i^aj, |Vxil<-, 
j 

where C is a universal constant. We define XR,j{^) = Xji^/^)- Then we obtain a new partition 
of unity XR,j such that supp XR,j C JCrj , with 

ICrj = {Rx : X G Kaj} . 
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Let J = {j = (ji,j2) G : 1 < Ji,i2 < n} and KnR = {-nR,nR) x {-nR,nR). Then the 
family {K,pi^j)j^j is a covering of KnR, and is formed exactly of ii? squares. 

We restrict the partition of unity {xR,j) to the set KnR = {—nR, nR) x (—nR, nR). Let UnR be 
a minimizer of G^j^k^r o^^^ the space E^/j, i.e. Gi,^K„R{'^nR) = 'mp{b,nR). We have the following 
decomposition formula (using the pointwise inequality '^Xrj ^ X^X/jj = !)> 

Gb,K^j,{UnR) > (Gb,K„^iXR,jUnR) - \\ |Vxrj| ^^nij|li2(^„^)) • (2.25) 

We reformulate ()2.25p as follows, 

n I \ ^ Gb,Kr,R{XR,jUnR) , a(2 + a) , ,, , , ,,2 \ 

Gb,K„R[UnR) > 2^ I (^l^ay (l + Q)2 '^^-^nfl(Xfij'»nflj " || I VXi?j I ^^nij || l2(^^^) I . 

(2.26) 

Notice that each XR,j has support in a square of side length (1 + a)R, hence it can be extended 
to a function in the space -E'(i+a)fi introduced in (j2.20p . Therefore, using magnetic translation 
invar iance, we get, 

Gb,K„R{XR,jUnR) > mp{b, (1 + a)R) . 

We insert this lower bound into ()2.26p . To estimate the localization error we use that |Vx_R,j| < 
C/{aR). The support of |VxR,j| is contained in an (ai?)-neighborhood of a square of sidelength 
(1 + a)R and therefore, its area is of the order aR^. So we get from (j2.26p the following lower 
bound, 

T-.\ 2"^p(^' (1 + '^)^) 2^*11 ||2 '^(2 + «) ^ 

mp{b,nR)>n . n -\\uni\\L^(^K^^) + . Gb,K,,R{XR,jUnR) ■ (2.27) 

Here the constant G is independent of a and R. Since XR,jUnR has compact support in M^, and 
the lowest eigenvalue of the Schrodinger operator with constant unit magnetic field in L^(M^) is 
equal to 1, we get by the variational min-max principle that, 

Gb,K„RiXR,jUnR) > (6 - 1) / IXRdUnRl"^ dx > -(1 - b)R^\\UnR\\loo(^K^^^) ■ 

J KnR 

Inserting this into (|2.27p then dividing both sides of the resulting inequality by n^R?, we get 
(using a < 1) 



mp{b,nR) mp{b,{l + a)R) / C 



_ ^ ■■■^^-^■^^^■^ _ + 3<.(1 - .)) ||u„„|li«,K-,„, . (2.28) 



We select eg such that < eq < l/yl + C^jg^, where Cmax is the universal constant from 
Proposition 12.61 In this way, we get that ||n„,R||2,°o(^^^) < CmaxVl — b for all 6 G (1 — eo, 1). 
We infer from ([Qop . 

rupib, nR) ^ rnpib,il + a)R) _ ^ f 4, + a(l - 6) Vl - 6) , (2.29) 



n2ii2 - (l + a)2i?2 \^a^2 

for ah 6 G (1 - eo, 1), > 2, n G N and a G (0, 1). Theorem [2T] and Proposition [2l\ together tell 
us that the function mp{b,£)/[i'^] has limit g{b) as £ — )• 00. Therefore, letting n — )• 00 on both 
sides of (|2.29p gives us. 
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We select a = 'jii^^jTjT^- It results from ()2.30p that 



9ib)> 



i?2 R 



for all b E (eo, 1), R> 2/-^/eo and some universal constant C. This finishes the proof of Proposi- 
tion [231 □ 



2.1.3. The Abrikosov energy. 

Recah the notation that if i? > 0, = {—R/2, R/2) x {—R/2, R/2) is a square of side length 
R. In this section, we assume the quantization condition that \Kr\/{2tt) is an integer, i.e. there 
exists N £'N such that, 

R^ = 2ttN. (2.31) 
Recall the definition of the space in (|2.20p . We denote by the operator, 

PR = -{V-iAof in L\Kn), (2.32) 

with form domain the space Er introduced in (j2.20p . More precisely, Pr is the self-adjoint 
realization associated with the closed quadratic form 

EB Bf^ QrU) = ||(V - iAo)/||i.(^^) . (2.33) 

The operator Pr being with compact resolvent, let us denote by {/ij(PR)}j>i the increasing 
sequence of its distinct eigenvalues (i.e. without counting multiplicity). 

The following proposition may be classical in the spectral theory of Schrodinger operators, 
but we refer to [2] or for a simple proof. 

Proposition 2.9. Assuming R is such that \Kr\ E 27rN, then the operator Pr enjoys the fol- 
lowing spectral properties: 

(1) fiiiPR) = 1 andfi2iPR) > 3. 

(2) The space Lr = Ker(PR — 1) is finite dimensional and dimL^ = \Kr\/{2-k) . 

Consequently, denoting by Hi the orthogonal projection on the space Lr (in L'^{Kr)), and by 
112 = Id — III, then for all f E D[Pr), we have, 

The next lemma is a consequence of the existence of a spectral gap between the first two 
eigenvalues of Pr. It is proved in [161 Lemma 2.8]. 

Lemma 2.10. Given p > 2, there exists a constant Cp > such that, for any 7 E (0, ^), > 1 
with \Kr\ E 27rN, and f E D{Pr) satisfying 

Q/?(/)-(l + 7)ll/lli2(^^)<0, (2.34) 

the following estimate holds, 

\\f-n^f\\L.iKR)<CpVl\\f\\LHKR)■ (2.35) 
Here Hi is the projection on the space Lr. 

We introduce the following energy functional (the Abrikosov energy), 

Fr{v) = (^-\v\^ + ^\v\''^ dx. (2.36) 
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The energy Fji will be minimized on the space L/j, the eigenspace of the first eigenvalue of the 
periodic operator Pr, 

Lr = {u £ Er : Pru = u} 

= {u£Er : {d^,+id^^+]^{xi+ix2)^u = Qi} . (2.37) 

We set, 

c{R) = m.{Yi{FR{u) : u € Lr} , (2.38) 

In the next theorem, we exhibit a connection between the ground state energies mp(b, R) and 
c{R) . In the regime where the parameter b is close to 1 , Theorem 12.111 is an improvement of 
Proposition 4.2 in [2]. 

Theorem 2.11. Let mp(b,R) and c{R) be as introduced in (|2.21|) and (|2.38p respectively. For 
all b > and R > 0, we have, 

mp{b,R) < [l-b]lc{R). 

Furthermore, there exist universal constants cq £ (0, 1) and C > such that, if R > 2, 6 > 1 — eo, 
and < cr < 1/2, then, 

mp{b, R)>[1- b]l ^(1 + 2a)c{R) - Ca~^{l - bfR"^^ . 

Proof. We start by proving the upper bound on mp{b, R). Let v G Lr be a minimizer of Fr, i.e. 

1 /2 

Fr{v) = c{R). Let / = [1 — b]_l v. Then / G Er and consequently, 

mp{b,R)<GKnU)- 
Using that v is an eigenfunction for the operator Pr, a simple computation gives, 

GK^f) = [1 - b]+ (-(1 - 6)b|2 + 1[1 - b] + \v\^^ dx 

= [1 - b]lc{R) . 

We prove the lower bound on mp{b,R). Let u G Er be a minimizer of Gkr, i-e. Gxniu) = 
mp{b, R). Since we do not know that u G Lr, we approximate u by its projection on the space 
Lr. Actually, we infer from mp{b,R) < the following inequality, 

/ (|(V - iAo)u\^ - \uf) dx<0. 

JKr 

Let 7 = [i - 1]+. We infer fr om Lemma [2.101 that . 

pin - u\\l4^k^) < G^\\u\\L2(^Kn) , 
where Hi is the projector on the space Lr. We deduce from this inequality the following estimate, 

\\u\\l4(Kr) > W^iuWlHKr) - C^/^\\u\\L^Ka) ■ 
This gives us, for some new constant C and for all a G (0, 1), 

Ml.^Kn) > (1 - ^)l|nin||i4(^^) - Ca-Vhllt2(^,) . (2.39) 

Now, we use the bound on ||n||j;^K)(iR2) given in Proposition 12.61 By selecting eo such that 
Cmax Veo/ \/l - eo < 1, we get for ah 6 > 1 - eo, 

M\.^Kn) ^ (1 - ^)l|nin||i4(^^) - Ca-^^^R^ . (2.40) 
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Using Proposition 12.91 and variational min-max principle, we also get, 

/ (6|(V -iAo)n|2 - l-up) > (6- 1) / \Uiu\'^ dx . 

JKr JKr 

This estimate, togther with that in ()2.40p give us the following lower bound, 

GKn{u) > j (^-(1 - b)|ni'u|2 + ^(1 - fT)|ni'u|^^ dx - CaS^R'^ . 

Recall that u is a minimizer of Gkh, and that Proposition 12.61 tells us that u = if 6 > 1. 
Therefore, the aforementioned lower bound is the same as, 

mp{b,R)>j^ (^-[l-b]+\Uiu\^ + ^{l-a)\Uiu\^^ dx-Ca-^j^R^ . (2.41) 

We introduce a function v £ Lr such that, 

[1-6]+^. 



1-a 

Notice that v = if b > 1, since u = in this case. Notice that v is constructed so that the right 
hand side on (|2.4ip becomes equal to, 

i — (7 

Since v G L/j, then Fpi^v) > c{R). Also, since < a < 1/2, then (1 - a)^^ < 1 + 2a. 
Consequently, we get, 

mp{b,R) > {1 + 2a)[l - b]lco{R) - Ca-^-f^R^ . 

Recalling the definition of 7 = — 1] + , the last estimate is nothing but the estimate of Theo- 
rem EJH □ 

As consequence of Theorems I2.1l and l2.11l and Proposition l2.8| we prove the following theorem. 

Theorem 2.12. Let the functions g{b) and c{R) be as defined in (|2.5p and ()2.38p respectively. 
The limits 

lim , and lim , \ 

b-s>l_ (1 - 6)^ R->oo Ar 

exist and are equal to a universal constant E2 G [—^,0). That is, 

E, = hm ^ = lim . (2.42) 

R^oo \Kr\ 6^1_ (1 - 6 2 ^ > 

Remark 2.13. The result of Theorem 12.121 appears first in [21 Theorems 1 & 2]. Our proof is 
based on a different approach from j^. Also in [21 p. 200], the proof of the inequality, 

lim sup 7-^—7 < liminf ■ ^^^^ 



R^J- \Kr\ - b^l. (1-6)2' 

\KR\/{2-K)m 

seems to rely on an estimate of the type we give in Proposition 12.81 which is missing in [2]. 
Proof of Theorem \2.12[ Using Theorems 12.71 and I2.11| we may write for all b G (^^ 1) and a G 

and 
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Here mo{b,R) is introduced in (|2.4p . We get by dividing both sides in (|2.43p by (1 — 6)^ then 
taking hminf, 

'^^^ <liminf4g^. (2.45) 



(1 - 6)2 - R^oo i22 • 

Consequently, by taking hmsup on both sides of (|2.45p . we get, 

b-5>l_ 

hm sup < Hm inf ^ . (2 .46) 

(1 - 5)2 R^oo B? 

Using the upper bound for mp{b,R) in Proposition 12.81 infer from (|2.44p . 

c(R) ^ 1 { g{b) , C , ,,,,2 



Let e > be given. By definition of Hmsup, there exists 6 > such that, if 1 — (5 < 6 < 1, 

gjb) . gib) 

n - limsup- 

(1 - by b^i_ (1 - b) 

b=l- we get for all R > (5~'^/^ 



then —K < limsupT —k + £. Inserting this upper bound into ()2.47p . then selecting 



<- (>fj- IT^ - - - ■ 

Taking successively limsup, then lim and lim on both sides of the above inequality, we get, 

c{R) gjb) 
hm sup —5- < hm sup —r . 

Combining this inequality with that in (|2.46p . we deduce that 

E2 := hm = lim sup -/^^ • (2-48) 

We return to (|2.47p and select R = {1 — 6)~^/^, so that — )■ 00 as 6 — )• 1_. That way we get 
(for this choice of i? = R{b)) 

Since c{R)/R? has a limit as i? — t- 00 which is given in (j2.48p . we get by taking lim inf on both 

b— 5>1_ 

sides of (IZiOjl . 

y gjb) r c{R) . , gjb) 
hm sup --77 = hm < hm mf --7 . 

(1 - 6)2 R^oc i?2 - (1 _ 5)2 

It results from the estimate ()2.6p in Theorem 12.11 that E2 G [— a]. The constant a is 
universal and a € (0, \). This finishes the proof of the theorem. □ 

2.2. Three-dimensional limiting energy. Let us again consider a constant 6 > 0. For any 
domain D C M'^, we define the following Ginzburg-Landau energy, 

Fl^{u) = j (^6|(V-iF)n|2- |n|2 + l|u|4^ dx. (2.50) 

We will sometimes omit the parameter 6 from the notation and write instead of F^g. 

Here F is the canonical magnetic potential, 

1 „ 
F(a;i,a;2,a;3) = -(-X2,xi,0) , V x = (xi, X2, X3) G M . (2.51) 
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We introduce the ground state energy, 

Mo{b,R)= inf Fj^^^iu). (2.52) 

Surprisingly, we find that the thermodynamic hmit of the functional F^ is equal to the 
corresponding two-dimensional limit. 

Theorem 2.14. Suppose R > 0, Qr = {-R/2, R/2) x Kr C R-^ and Kr = {-R/2, R/2) x 
{—R/2, R/2) C M^. The following statements hold true. 

(1) Ifb>l, then for all R>0, Mo{b, R) = . 

(2) There exists a universal constant M > such that, for all b > and R > 0, we have, 

Rmo{b,R) < Mo{b,R) < {R - 2)mo{b,R) + M , 



where mo{b,R) is the ground state energy introduced in ([2? 
(3) 7/0 < 6 < 1, then 

fl^oo \Qr\ 



where g{b) is the constant from 

Proof. We start by proving the statement corresponding to 6 > 1. Using the configuration u = 
as a test function, we get obviously that 

inf FS°(n)<0. 

Let u G Hq{Qr). To get a lower bound, it is sufficient to prove that Fq^{u) > 0. We extend u 
to a function u e i7^(M^) by setting 5 = in M'^ \ Qr. Then, 

i^a^(") = ^3(&l(V-^F)S|2-|S|2 + l|5|4^ dx. (2.53) 

Recall that the bottom of the spectrum of the magnetic Schrodinger operator (—V — zF)^ in 
L^(M^) is equal to 1. Using the variational min-max principle, this gives the following lower 
bound, 



/ \{V - iF)u\^ dx > [ 



Inserting this into ()2.53p . we get that F|°(n) > 0. 

We prove the statement corresponding to < 6 < 1. First we mention that if x = {xi,X2,X3) £ 
R^, then we write x = {x±,X3) and Va^^ = {8x1,8x2)- Let u S Hq{Qr). Notice that 



|(V - iF)u\' = \{Vx^ - iAo)u\' + \dx,u\' > \{Vx^ - iAo)u\ 
This gives the following lower bound. 



fR/2 

F^l{u)> / GKM-^^3))dxs, 

J -R/2 



where Kr = {—R/2, R/2) x {—R/2, R/2) and is the functional introduced in ()2.ip . Since 

u S Hq{Qr), then u{-,X3) E Hq{Kr), for X3 S {—R/2, R/2) almost everywhere. So we can write 
GKR{u{-,X'i)) > mo{b,R). 

Consequently, we deduce that 

Mo{b,R) > RmQ{b,R) . 
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To get a matching upper bound, we consider the test function f{xi,X2,xs) = X2)X-r(^3)) 
where Uh{xi,X2) is a minimizer of Gkji and the function Xni^s) satisfies, 

Xr£C^{R), 0<Xi?,<l, suppxRG [-i?/2,ii/2], ;^^ = 1 in + 

and \x'r\ ^ M in M, for some universal constant M. 

We estimate the energy FQ^{f). We get, using the pointwise bound x% ^ X/ji 

if/2 ^ /■_R/2 ^ 



F^°(/)<f/ \u,\^dxidx^f \x'R{xz)?dx-i + GKM 

yJKa / J-R/2 



Xr{x3)\ dx3. 

R/2 

By (|2.3p we have \uh\ < 1. Also, using the properties of XR g^t, 

/ \x'R{xs)\^dx3<2M, R-2< \xR{xs)\Ux-i<R. 

J-R/2 JR 

Thus, we get, 

F3D(/)<2M + (i?-2)Gx,M. 

Consequently, we obtain, 

Mo{b, R) <2M + {R- 2)mo{b, R) . 

This proves statement (2) in Theorem 12.141 

The last statement in Theorem 12.141 results straightforwardly from the inequality in the second 
statement. Actually, we divide both sides of the inequlity by R^ then we take i? — )• oo. □ 

3. A PRIORI ESTIMATES 

The aim of this section is to give a priori estimates on the solutions of the Ginzburg-Landau 
equations ()1.3p . Those estimates play an essential role in controlling the error resulting from 
various approximations. 

The starting point is the following L°°-bound resulting from the maximum principle. Actually, 
if (■0, A) G H^{n; C) X H^.^^iM.^) is a solution of ([O]), then 



L°°(n) 



<1. (3.1) 



The set of estimates below is proved in |15| Theorem 3.3 and Eq. (3.35)] (see also |22) for an 
earlier version). 

Theorem 3.1. 

(1) Suppose 1 < p < 6. There exists a constant Cp > such that, if k > 0, H > and 
(-0, A) G H^{n) X i7i;^ p(M3) is a solution of ([TS]), then 

I _|_ _|_ ^2 

||A - F\\w2.p(^n) < Cp — M\L^in)U\\L^in) ■ (3.2) 

(2) There exists a constant C > such that, if k > 0, H > and {tp, A) G H^{i})xH\^^-p{M.^) 
is a solution of (jl.3p . then 

||curl(A-F)||i2(^) < ^\mL^in)\mmn)- (3.3) 

Theorem 13. II is needed in order to obtain the improved a priori estimates of the next theorem. 
Similar estimates are given in [22) . 
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Theorem 3.2. Suppose that < Amin < Amax- There exist constants kq > 1 and Ci > such 
that, if 

and (V', A) G i?^($7;C) x if]j^ p(M3) is a solution of ([O]), then 

II (V - iKFA)^llp(^) < Ci V^||V^||i^(f,) , (3.4) 

||A - F||ty2.6(n) < Ci ^11 curl(A - F)||i2(ffi3) + -^||V'||L6(t^)||V'|U-(Q) j , (3.5) 
||A - F||^M/2(n) < C\ (ll curl(A - F)||i2(M3) + -^||V^||i6(o)||^||L^(n)) ■ (3.6) 

Proof. 

Proof of dSai): 

Suppose the estimate (|3.4p were false. Then, there exists a sequence of points {Pn) C 0, 
two real sequences (/?«), (Hn), a sequence of solutions {^l'n,A.n) of (jl.Sp and a constant A £ 

[Amin, Amax] SUch that, 

|(V - iKnHnAn)tpniPn) 



OO, (3.7) 



Kn — )• +00 , — > A as n — )• OO . 



Let (S'ti 

= IIV'n||L°°(n)- Define the re-scaled functions. 



and 



V VKnHnJ 



A„(P„ + y/^/KnHn) - An{y) 

Aniy) = 



1/ sf^^nHn 

Two cases may occur: 



Case i.' KnHn dist(P„, 90) is an unbounded sequence. 
Case 2: \/KnHn dist(P„, Sfi) is a bounded sequence. 

In both Cases 1 and 2, we apply a compactness argument by using elliptic estimates then a 
diagonal sequence argument to select a limiting function. Actually, we refer to \15\ Section 4.3] for 
the detailed proof in the two dimensional case, and to \12\ Lemma 12.5.4] for a precise statement 
in the three dimensional case. 

Therefore, in Case 1, we can extract a subsequence, still denoted cpn, a function £ 
and a vector field F G C(]R^) such that 

tpn^if in C^{K), An^F in C{K), 
for any compact set ii' C M^. In particular, we get that, 

liV - iK.nHnAn)lpn{Pn)\ , / m 

^^i" , "^"^ ^ |(V - zF)(^(0)| , 

yl^'nJ^n \\iPn\\L°°{n) 

which is contradictory with p.7p . 

Similarly, in Case 2, we can extract a subsequence, still denoted a function ip £ 
and a vector field A G C(]Rij_) such that 

ifn^ip in C\K), An^A in C7(i^) , 
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for any compact set K C M'^ . In particular, we get that, 

i'^-!!;;-yf"'^'''U |(v-iA)^(o)i. 

which is again contradictory with (|3.7p . Therefore, the estimate (|ci.4p should be true. 
Proof of (I33]l : 

Let a = A — F. Since diva = 0, we get by regularity of the curl-div system see e.g. (or |12| 
Theorem D.3.1] for a statement of the result), 

l|a||i,6(R3) < C\\ curl all i2(iR3) . (3.8) 
The second equation in (|1.3p reads as follows, 

-Aa = Im(V^ (V - iKH A)i;)lQ . 
By elliptic estimates (see e.g. or \12\ Theorem E.4.2]), 

II«I|W'2,6(Q) < C(||a||2,6(K3) + ||Aa||^6(K3)) . 

The estimates in p.4p and (|3.8p now give, 



/ C 

||a||w'2,6(Q) < C II curla||i2(]R3) + \\'4'\\L<i{n)ML^{n) 
\ V kH 

Proof of IK^ : 

This is a consequence of the Sobolev embedding of W'^^^{n) into C^''^/^{n) and ([33]). □ 

The next result is a rather weak L°°-bound valid for all critical points of provided that 
the magnetic field strength H is close to Hc2- 

Theorem 3.3. Suppose that gi : — )• M4. is a function such that 

( \ . 5'i('^) , n 
51 (kJ — )• 00 , )■ U as K — )■ 00 . 

There exists a function §2 ■ — )■ M+ satisfying 

92(1^) — 00 as K — 7- 00 , 

such that if n> 1, \H — k\ < gi{n) and G H^{Q,;C) x ij|^jyp(M^) is a critical point of the 

energy in (jl.ip . then, 

II^IIl'-K) < 92{k) . 

Here 

uJk = {x€^1 : dist(x, (9f]) > (7i(k)/k} . 



Proof. The proof is very close to Theorem 2.1 in |13| . We give the details for the reader's 
convenience. 

Suppose that the conclusion of Theorem 13.31 were false. Then, we may find a number N > 1 
and sequences Hn, {tpn,-^n) such that, 

Kn — )■ +00 , — 1 as n — +00 , 

iV"' < ||V'n||L-(^.„) < 1, (3.9) 

and {ipn,-^n) is a solution of (|1.3p for {k,H) = {Kn,Hn)- 
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It results from Theorem 13.11 that ||A„ — F||^2,6(q) < Ci for some constant Ci > 0. Using 
the compactness of the embedding of W^'^{^1) in C^'"(ri) for < a < ^, we may select a 
susbsequence, still denoted and a function A G C^'°'{Q) such that 

A„^A in C^'^p. (3.10) 

It results from the estimate in p.3p and the assumption Hn — s- oo that 

curlA = l inO. (3.11) 

Let Pn G uj^^ be a point satisfying |'0„(P„)| = ||V'n||L°o(i.jK„)- Using ()3.9p . we deduce that, 

N^^ < \lP{Pn)\ < 1. 

We may select a new subsequence, still denoted {Pn;ipm -^n), and a point P G 17 such that 
Pn ^ P as n ^ +00. 

We define the following re-scaled functions. 



My) = . , , — rj- 



Since gi{nn) — ^ +oo as n — )• +oo by assumption, it follows that for any > 0, we may select no 
sufficiently large such that and are defined in {y G : |y| < R] for all n > uq. 

Notice that a„ is constructed so that diva„ = 0. We then infer from the equation of ipn the 
following equation, 

- A(pn - 2ia.n ■ Vifn + |a„p99„ = -^(1 - \ipn\^)'^n ■ (3.12) 

Consider R > 1. Using the definition of a^ together with ()3.10p and (|3.1ip we deduce that, 

a„^F(y) in C0'"(5(0,i?)), 

where F(y) = DA{y) satisfies curlF = 1. 

Since \ipn\ < 1 and |a„| < Cr in -6(0, 2R), we get by elliptic regularity that the sequence (</?„) is 
bounded in W^'^{B{0, R)). The compactness of the embedding W'^'^{B{0, R)) Ci'°(5(0, R)), 
a < ^, and a standard diagonal sequence argument give the existence of a function (p G C[^'°(M'^) 
such that 

ipn^^ in Ci'"(K), 
for each compact set C M^. Furthermore, the function ip satisfies, 

N^^ < ||v'||l^(r3) < 1 , 

and 

_(V-iF)V = (l-j¥'|')</'. 
Since curlF = 1, we may find a function (j) such that F = F + where F = (— 2:2/2, xi/2, 0). 



Setting u = e ^'''tp, we get that u G L°°(M ) and satisfies the equation, 

-{V -iFfu = {l-\u\^)u. (3.13) 

At the same time u satisfies, 

||n||ioo(K3) > iV-i . (3.14) 

But Proposition 12.5.1 in |12) tells us that the only bounded solution of ()3.13p is u = 0, thereby 
contradicting ()3.14p . Therefore, Theorem 13.31 is true. □ 
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4. Energy estimates in small cubes 

In this section, the notation Qi stands for a cube in of side length i > and whose axis is 
parallel to (3 = (0, 0, 1), i.e. Qi is of the form, 

Qi = [-112 + ai, ai + ^/2) x (-^/2 + 02, 02 + ^/2) x (-^/2 + 03, 03 + ijl) , 

where a = (ai, 02, 03) € M'^. 

If (i/;, A) G ^1(0; C) x Hl^^^iM.^), we denote by e{ij, A) = | (V - KHA)i;\'^ - K^jV'P + ^IV-jl 
Furthermore, we define the Ginzburg-Landau energy of (-0, A) in a domain 2? C ^2 as follows, 

£{ij,A;V)= [ e{t^,A)dx + {KHf [ | curl(A - F)|2 . (4.1) 

We also introduce the functional, 

£o{u,A;V) = J (^\{V - iA)u\^ - K^\u\'^ + dx . (4.2) 

If 2? = J7, we sometimes omit the dependence on the domain and write £o{ip, A) for Soiijj, A; Q). 

We start with a lemma that will be useful in the proof of Proposition 14.21 below and also in 
Section [6] of the paper. 

Lemma 4.1. Let < Amin < Amax- There exist positive constants C and kq such that if 

K 

and if (ipjA) G H^{Q;C) x i?^-^p(M'^) is a critical point of (jl.ip . f is a bounded function 

satisfying \\f\\L°°[n) < 1, supp/ C Qi, ftp € H^{Qi>) and Qi <Z Q is a cube of side length 
i G (0, 1), then the estimate, 

£o{fi^,A-Qi)>il-5)£oie-''^°fi^,F;Qi)-C{6K^+5-H^X) [ dx , 

jQi 

holds true for all 6 G (0, 1), k > kq and some real-valued function (f)Q G H^{^). Here 
A = {^Hf curl(A - F)||i.(j,3) + ^M^in)) ■ 

Proof. Let xq be the center of the cube Qi. Without loss of generality, we may assume that 
xq = 0. In this way, we reduce to the case where 

Qi^ = (-4/2,4/2) X (-4/2,4/2) X (-4/2,4/2) cn. 

Let (j)o{x) = (A(xo) — F(rEo)) -x, where F is the magnetic potential introduced in (|2.5ip . Invoking 
the estimate in (|3.6p . we get, 

\A{x)-V(j)o-Fix)\<C^i WxGQi. 

Let u = e-'^ofijj and A = A - V(/>o. Then Soifil^, A; Qi) = £q{u, A; Qi) and 

|A(x)-F(x)|<^^, VxgQ,. (4.3) 

We estimate the energy £q{u, A; Qi) from below. We start by estimating the kinetic energy from 
below as follows, 

|(V - iKHA)u\^ > (1 - 5)\{V - iKH¥)u\'^ + (1 - 5-^){kH)'^\(a - Y)u\^ , 
for any 6 G (0, 1). 
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Using the estimate in (|4.3p together with the assumptions -ff ~ k, ||/||oo ^ 1 we deduce the 
lower bound, 

Soifi^, A; Qe) > (1 - 6)£oiu, F; Q,) - C{6k^ + 6"H^\) [ dx . (4.4) 

□ 

Proposition 4.2. Let < Amin < Amax- There exist positive constants C,Rq and kq such that 
the following is true. 

Let £ > satisfy Rq < i < 1/2 and let 

H 

K 

Then, if {tp, A) G H^{Q; C) x ij^^^^ p(M'^) is a critical point of (|l.ip and Qi C Q is a cube of side 
length I, then, 



rSoii^, A; Q,) > g{H/K)K^ -c[1 + —]k 



2 



/or all K > 1. i/ere is the function introduced in (|2.5|) . 

Proof. We may reduce to the case where 

= {-1/2,1/2) X {-1/2,1/2) X {-1/2,1/2) cQ. 

Recall the parameter A introduced in Lemma [4.11 Using the bound Hi/'lloo ^ 1 and the estimate 
in ()3.3p . we get A < Ck^. We apply Lemma |4. II with / the characteristic function of Qi. After 
possibly performing a gauge transformation, we may assume the following lower bound, 

£o{i;,A;Q,) > {1 - 6)£o{i^,F;Qe) - C{6 + 6-H^)K^mml^(^Q^^ , (4.5) 

where C > is a constant only depending on Q. 

We estimate the energy i5o('0) F; Qf) from below. Let b = H/k and R = I^/kH . Define the 
rescaled function, 

ip{x) = iP{x/VkH) , y X £ Qr . 

Recall the functional F'q^ introduced in (I2.50p above. It is easy to check that, 
£,{^,F;Qe) = -^l^ (l(V - .F)^!^ - -^H^ + ^H^) dx 

^ -.F'^M- (4.6) 



b^/KH 

We still need to estimate from below the reduced energy F|°('i/'). The probl em is that if) is not 
in Hq{Qji). So, we introduce a cut-off function XR ^ C^(M'^) such that 

< XiJ < 1 in , suppxR cQr,, Xij = 1 in . 

In particular, the function XR can be selected such that |Vx_r| < C for some universal constant 
C. Let u = XR^^- It is easy to check that, 

F3D(V^) > / (b\xR{V - iF)V^|2 - IxfiV^p + IxRi^lA dx+ [ {xl - i)|V^P dx 
JQr ^ ' JQr 

>^Q°(")+ / {\^XR?+XR-m?dx-2\{{V -iY){xR^J),{VxR)i^)\■ 
JQR 
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The estimate in ()3.4p tells us that |(V — iF)'ip\ is bounded independently of k and H. Also, by 
using the bounds \u\ < 1, |Vx_r| < C and the assumption on the support of XR^ it is easy to 
check that, 

After recalling the definition of Mo(6, R) introduced in ()2.52p . we get, 

> Mo{b,R) - CR^ ||Vllioc(Q,) . (4.7) 
We get by collecting the estimates in (I4.5|) - (|4.7p that, 

^ £o{^,A;Qe)>{l-6)—^Mo{b,R)-r{K), (4.8) 



\Qi\ 'u^V^ 

where 

r{K) = c({5 + b-H^y^ + — ^ ) II Vllioo.o . . (4.9) 



R^ 

Theorems O and EH together tell us that Mo{b, R) > R^gib) for all b G [0, 1] and R sufficiently 
large. Here g{b) is introduced in (|2.5p . Therefore, we get from ()4.8p the estimate, 

^ fo(V',A;Q,) > (1-5)— ^5(&)-r(Ac). (4.10) 



IQ^I '-^^ - ^ ^6^3^ 

We choose 5 = £ so that 

r(K) = 0((^ + (^k)-1)k2). 
After recalling that R = i\J kH and b = H/ k, we finish the proof of the proposition. □ 

Proposition 4.3. Let < Amin < ^max- There exist positive constants C,Rq and kq such that 
the following is true. Let £ > satisfy Rqk~^ 1^ £ 1^ 1/2 and let 

H 

Then, if [ip^A) € H^{Q;C) x i/^j^p(]R'^) is a critical point of p.ip . and d Q is a cube of 
side length £, then, 

^ -M^, A; Q,) < g{H/K)K^ + c(£ + ^\k\ 



Here g{-) is the function introduced in ()2.5p . and £q is the functional in ()4.2 



Proof. After performing a translation, we may assume that, 

= {-£/2,£/2) X {-£/2,£/2) x {-£/2,£/2) cfl. 

As explained earlier in the proof of Lemma 14. 1| we may suppose, after performing a gauge 
transformation, that the magnetic potential A satisfies, 

C£ 

|A(x) -F(x)| < VxGQf, (4.11) 

where F is the magnetic potential introduced in ()2.5ip . 

Let b = H/k, R = £\fKH and ur € HI{Qr) be a mmimizer of the functional F^D introduced 
in (I230D . i.e. F^^{ur) = Mo{b,R) where Mo{b,R) is introduced in ([232]) . 

Let xr £ C^(M^) be a cut-off function such that, 

< Xfl < 1 in , supp XR C Qr+i , Xfl = 1 in . 

and |Vx_r| < C for some universal constant C. Let riR{x) = 1 — Xr{xV kH) for all x G M'^. We 
introduce the function, 

(p{x) = 1q^{x)ur{xVkH) + ■qR{x)ilj{x) , V x G r2 . (4.12) 
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Notice that by construction, if = in Q\ Qp , i . We will prove that, for all 6 G (0, 1), 

£{ip, A; Q) < f (V, A; Q\Q,) + {1 + S)-^Mo{b, R) + ro{K) , (4.13) 

b\' kH 

where Mo(6, -R) is defined by ()2.52p . and for some constant C, ro(K) is given as follows, 

ro(^,) = Ci^ (^5 + + j-^ k\ (4.14) 

Before proving (j4.13p . we explain how we get the result of the proposition. Actually, by 
definition of the minimizer (■0, A), we have, 

^(V',A) <£{ip,A-n). 

Since A; il.) = A;i} \ Qi) + <fo(V') A; Qe), the estimate ()4.13p gives us, 

£o{i;,A;Qe) < {I + 6)-^Mo{b, R) + Mk) . 
by nH 

Dividing both sides by \Q(\ = £^ and remembering the definition of ro(K), we get, 

^ r£{^, A, Qt) < ^\'^lL MQ{b, R) + c(6 + 6^^f + ^] k\ (4.15) 



\Qi\ " ' - bi^^/l^ 

Theorems O and EH together tell us that Mo{b,R) < R^g{b) + CR^ for aU b G [0,1] and R 
sufficiently large. We substitute this into (|4.15p and we select 5 = I, so that ro(K) = K^O{l + 
[(.k)~^^. Recalling that R = i\J kH , b = H/k, we get, 

t)-M^, a, Q,) < g{b)K^ + K^O{£ + {Ik)~^) . 
\Qe\ 

This establishes the result of Proposition 14.31 
Proof of (ITOD : 

Recall the Ginzburg-Landau energy £q defined in ()4.2p . We write, 

£o{^,A;n)=£i+£2, (4.16) 

where 

£i=£o{v^,A;n\Qe), £2 = £o{v^, A;Qi) . (4.17) 
We estimate each of £1 and £2 from above. Starting with £1, we write, 

£i=£oi^l^,A;n\Qi)+nii;,A), (4.18) 

where 

n{^P, A)= [ ({tjI - 1) (|(V - iKHA)^P\' - ^\^P\') + |^Vr?«|2 + ^(^ - 1)1^^14 



+ 23f?(?7R(V - iKHA)i), ipVr]R) dx . 



Using that ijji = 1 in \ , 1 together with the estimates iVry/jj < CVkH, < 1 and (|3.4p . 

we get, 

\n{'ilj,A)\ <CfK. (4.19) 
Inserting ()4.19p in (j4.18p we get the following estimate, 

£1 <£o{ip,A;Q\Q(;) + Ce^K. (4.20) 
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We estimate the term £2 from ()4.17p . We start by observing that for any 6 G (0, 1), we have the 
following upper bound, 

^2 < ^ |(1 + <5)|(V - mOT)c^|2 - k2|(^|2 + y |<^|'| dx 

+ {1 + 5-^)k^H'^ f |A- F|2|(^|2dx. (4.21) 
JQi 

We use the estimates ()4.1ip and < 1 and we get, 

£2<{l + 5) j - iKH¥)^\^ - + y dx + C5-H\^ . (4.22) 

Since (p{x) = n/j(x\/ kH) in Qi and R = I^/kH, we get by performing a change of variables, 
(|(V-ZKi?F)v.|2-KVl' + Yl^r) dx 

^ F'^'iM, (4.23) 



where is the functional from (|2.50p (with b = H/k). Inserting (j4.23p into (|4.22p . we get 

1 



£2<il + d)^^F^^^{uB) + C5fK^ + C5-H^K^ . (4.24) 
by kH 

Inserting (|i:20]l and (filMl) into ()irT6]l . we deduce that, 

^o(^, A) < £:o(V', A; !^ \ Q^) + (1 + 5)^^F^° (ur) + CSI^k^ + C6-H^k^ + C^^^ . (4.25) 
Recalling the definition of £{ip,A; •) in ()4.ip . we see that ()4.25p is sufficient to finish the proof 

of Km . □ 

We conclude the section by giving an upper bound on Eg.st(/^, H). 

Proposition 4.4. Let < Amin < Amax o^iT'd let 5 > 0. There exist positive constants C and kq 
such that if 

H 

then the ground state energy Eg.st(K, ff) in (|1.4p satisfies, 

Eg.st(K, iT) < 5(F/k)|J7|k2 + Ck^+'^ . 

Proof. Let ii = i?(K) be a positive parameter such that l^i^^KasK— )-oo. We will choose 
ii as a power of n at the end of the proof. Recall the magnetic potential F introduced in (|2.5ip 
and the ground state energy Mo(6, R) in (j2.52p . 

Let (Qj) be the lattice generated by the cube Qfi/^/^- I^et I = {j : Qj C 0}. We define 

= Card (X). Then, as R/\/kH 0, N satisfies. 



"=(v3j 

Let b = H/k, and u a minimizer of the functional in ()2.50p . i.e. 



Mo{b,R)= [ (b\{V-iF)u\ 

Jqr V 



+ 2 1''^!^ 1 dx . 



v{x) 
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Here recall the notation that if r > 0, then Qr = (-r/2,r/2) x (-r/2,r/2) x (-r/2,r/2). 
Since u E H^(Qji), we extend u to all by 'magnetic periodicity' such that, 

u{xi + R, X2,X3) = e^^'^^'^u{xi,X2,X3), u{xi,X2 + R, X3) = e~^^'-^^^'^u{xi,X2,X3), 

u{xi,X2,X3 + R) = u{xi,X2,X3) . 

For all X G ri, define, 

u{x\/kH), X G Uj^xQj 
^0 if not 

Since the original u G Hq{Qr), this defines a w G H'^{Q). 

We will compute the energy of the configuration (f,F). By periodicity, we get, 

£{v,F)=Nx[ ( \{V - iKHF)v\^ - K^lvl"^ + —\v\A dx . (4.26) 

A change of variables gives us, 

/ f |(V - ^.HF)v\^ - + ^\vA dx = . 

Inserting this in ()4.26p . we get, 

We know from Theorems O and EH that Mo{b,R) < g{h)R^ + CR^ for ah h G [0,1] and R 
sufficiently large. Using this together with the estimate on the number N , we get. 

We select R = k~^^^ . Recalling that b = H/k, this choice of R and the aforementioned upper 
bound finishes the proof of the proposition. □ 



5. Proof of main theorems 

Proof of Theorem 

In light of Proposition I4.4| we only need to establish a lower bound for the ground state energy 
Eg,st(K,//) in (HaD. 

Let I = Then k"^ < 4 < 1 as K — )• 00. Consider a lattice {Qj/)j of M'^ generated by 

the cube, 

= {-£/2,£/2) X {-£/2,£/2) x {-£/2,£/2) . 

Let I={j : Qj/ C n}, iVs = CardX and Q^, = VL\{[Jj^xQj/)- Notice that = \9\£~''^ +0 {£~'^) 
as K — 00. 

If (V') A) is a minimizer of (jl.ip . we have, 

Eg,,t(K,F) =fo(^,A;0«) + £:o(^,A;0\J7«) + (K/7)2 /" | curl A - Fj^ dx , (5.1) 

where, for any D C 0, the energy So^tp, A; D) is introduced in (|4.2p . We may write, 

Eg.st('^, H) > fo(V', A; n,) + £:o(V, A; J7 \ J],) . (5.2) 
Notice that [fi \ fi^l ~ £\dQ\ as k — )■ 00. Thus, we get by using the estimate in (|3.4p . 

\£o{tp,A;n\n^)\ <C£k^ = Ck^/^ . (5.3) 
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To estimate £q{iP, A; Ok), we notice that, 

Soiil), A; Ok) = ^ £Q{ij, A; Qj- ^) . 

Using Propositions 14.21 we get, 

fo(V', A; Ok) > [g{h)K^ - Ck^'^) ^ = Mg{h)K^ - Ck^'^ . (5.4) 

Inserting (|5.4p and (|5.3p into (|5.ip we deduce that, 

^^A^^,H)>g{b)\n\K^-CK^'^, 

thereby finishing the proof of Theorem 11.11 □ 

Proof of Theorem \1.2l It suffices to prove Theorem II .21 in the case where i < In order to 

reduce to this case consider a cube Qi with sidelength i > We can subdivide this cube 

into {i/£')^ identical cubes of side-length £' G [k~^/^/2, k~^/^). If the result of Theorem 11.21 is 
valid for each of these smaller cubes, the statement for Qi follows. 

Therefore, let i be such that <^ i < k^^/^, and Qg a cube of side-length i with a side 
parallel to the external magnetic field /3. Let (V', A) be a solution of ()1.3p . Then ip satisfies, 

- (V - iKFA)^^; = (1 - |V'P)V' in ^- (5.5) 

We multiply both sides of the equation in ()5.5p by ■0 then we integrate over Qg. An integration 
by parts gives us, 

[ (|(V-iKi7A)V'P -kVP + ^VI^) /" 1^ ■ (V -inHA)i;lpda{x) = 0. 

JQi JdQt 

Using the estimates ()3.ip and (j3.4p . we get that the boundary term above is 0{k£'^). So, we 
rewrite the above equation as follows, 

-^K^/ \i;\Ux = £o{^P,A;Qe) + 0{Kf). (5.6) 
2 Jq, 

Using Proposition 14.21 and the assumption ^ i, we conclude that 

IV'I^ dx < -g{H/K)\Q,\K^ + c(^£+j^^ K^£\ (5.7) 

If (■0, A) is a minimizer of (jl.ip . then (j5.6p is still true. We apply in this case Proposition 14.31 
to write an upper bound for £o{ip, A; Qi^)- Consequently, we deduce that, 

^K^l \i;\Ux>-g{H/K)\Q,\K'-c(^£ + ^^KH^. (5.8) 



Combining the upper bound in (|5.8p with the lower bound in ()5.7p and using that by assumption 
£ < /t~^/^, finishes the proof of Theorem 11.21 □ 

6. Additional estimates 

The aim of this section is to give additional estimates on the energy which are particularly 
interesting in the case where the magnetic field satisfies 

H = K + o(k,) as K — )• oo . 

These estimates will be used in [18] to prove the asymptotic formula in (|1.6p . 
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6.1. A lower bound. 

We will prove Theorem 16. II below, whose statement requires some notation. Let D C he a. 
given open set such that there exists a subset D of having smooth boimdary and D = DPiil. 
For all a > 0, we assign to D the following subset of Q, 

Da = {xen : dist(a;, D) < a} . (6.1) 

Let g : ]R_|_ — )• ]R_(_ be a function such that g{K) — )• oo as k — )■ oo. We put 

u^K = € ^ ■ dist(x,0r2) > g{K)/K} . 

Theorem 6.1. Suppose that the magnetic field H is a function of k such that, 

H H 
1 < lim inf — < lim sup — < oo . 

Let K 3 R+ I—)- a{K.) E be a function satisfying lim a(K) = 0. Then, for any solution 

(■0, A) G H^{^};C) X iir^i^ p(]R^) of ([TT]) and any function h G C^{0,) satisfying \\h\\i^(^Q) < 1 
and supp/i C Da D lo^^, the following asymptotic lower hound holds, 

£o{hij,A) > E2\D\[k - H]l + o^max (k, [k — if]^) ^ , as k — ;> 00 . (6.2) 

Here E2 < is the universal constant introduced in ()2.42p . and £q is the functional introduced 
m (lOl. 



The proof of Theorem 16.11 will be split into several lemmas. One ingredient to control the 
errors resulting from various approximations is an L^-bound on the order parameter given in [5] 
(also it is obtained in |12| Chapter 12] by a different method). We state this bound below. 

Lemma 6.2. Suppose the magnetic field H is a function of k and satisfies, 

H , H 
1 < lim inf — < lim sup — < 00 . 

There exist positive constants C and kq such that, if A) G H^{^}; C) x i^T^j^ p(]R'^) is a solution 
of dO]), then. 



1 



At 



H \ 



1^1^ dx < C max 

JO. 

for all K> HQ. 

Proof. This is a combination of the conclusions of Lemma 3.2 and Theorem 3.3 in [3]. □ 

The next lemma is taken from \12\ Lemma 10.33], which, together with Lemma 16.21 give a 
good estimate of || curl(A — F)||j;^2( 



Lemma 6.3. There exists a constant C > such that, if {iJj,A) G H^{Q.;C) x ij^.^p(M^) is a 
solution of ()1.3p . then, 

II curl(A - F)||^2(K3) < ^U\\%(^n) > (6-3) 

for all K> Q and H > 0. 

Lemma 6.4. Suppose the assumptions in Theorem lg.il hold. There exist positive constants C 
and kq such that, i/(0, A) G H^{Q;C) x if^Jj^^ p(M'^) is a critical point of (jl.ip . f is a continuous 
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function satisfying \\f\\L°°{n) ^ 1; supp/ C Qc and Qi d Q, is a cube of side length I G (0,1), 
then the estimate, 



2 TT2e4 



£o{f'il^,A) > {1-6) — [k-H]1 [{l + 2a)c{NKH)-Ca-\K- HYH 



C{6k' + 5-H^\) I \fi^\^dx, (6.4) 



holds true for all 5 G (0, 1), a G (0, 1/2) and k > kq. Here, 



A = {kHY ( II curl(A - F)||i,(j,3) + 11^^1126(^7) 



(6.5) 



Proof. Using Lemma [4. H we get a real- valued function (f)Q such that, with u = e ^'^^ fij:, we may 
write, 

£o{f^,A.)>{l- 5)£q{u,¥) -C{5k^ + 5^^l^\) f \fi)\'^ dx . (6.6) 

JQi 

We estimate the energy £q{u,Y) from below. Notice that suppu C Qi- By defining R = 1\[kM 
and the re-scaled function, 

V X G Qb. , v(x) = u{x/\/ kH) , 
we get that v G Hq{Qji). An easy computation gives, 

£oiu,F) = ^l^ (l(V - ^F)v\' -if\vn dx . 

Recall the definition of MQ{b,R) in (|2.52p . With this notation, we get. 

The conclusions in Theorem 12.141 Theorem 12.111 and Proposition 12.71 altogether give. 



(6.7) 



Mn 



H 



,R] >R 



H 



(1 + 2a)c{R) - Ca-^ ( ^ " ~") ^ 



Recall that R = i\J kH . Substitution into (|6.7p gives, 

£q{u,¥) >-^[k- H]\ f(l + 2a)c{lV^) - Ca'^in - H)'^H^i^) . 
Insersting (|6.8p into (j6.6p finishes the proof of Lemma [6.41 



(6.8) 

□ 

Proof of Theorem \6.1[ Consider two parameters i and a in the interval (0, 1). The parameters i 
and a will be chosen as functions of k such that 

i ^ , a— )-0, asK— )-oo. 

We cover M'^ by cubes (Qi(xj^a))j, where for all j G Z^, y = {yi,y2-,yz) £ a-i^d £ > 0, we 
define, 

3 



cj> = (i-a)j, QKy) = n 



k=l 



Let {fj) be a partition of unity in M'^ such that. 



Vk- ^,yk + ^] ■ 



C 



^ /J = 1 , supp/j C <9i(2;j,a) , |V/j| < — , 



for some universal constant C. 
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Defining tlie re-scaled functions, 

we get a new partition of unity {fj/)j such that each fj^£ has support in a cube of side length I 
and, 

j 

Let J = {j : DaH supp/j/ ^ 0} and Me = CardJ". Then we know that 

/?3 



We have the localization formula. 



£o{hij,A) = y2£oif,,M,A) - V / |V4,|2|/iV^|2dx. (6.9) 



Using that h < 1, (|6.9p gives us the following lower bound 



Soihi;, A) > ^ SoihM, A) - / IVP dx . (6.10) 



Many of the error terms will be controlled by the following parameter, 

K, 



= max I K 



-1/2 



H \ 



(6,11) 



Recall the parameter A from Lemma |6.4[ Using Lemmas 16.21 and 16.31 we get the following upper 
bound on A, 

A < Ck^C"^/^ . 

We apply Lemma [6.41 with / = fj^£ to bound from below each term i5o(/j,£^^i A). This gives us. 



£o{hi;, A) > Mi{l -6) — [k- H]\ (1 + 2a)c{iV^) - Ca~^{K - HfH^i^ 



\^p\'^dx, (6.12) 



for ah 6 G (0,1) and a G (0,1/2). 

Theorem I2.12l tells us that if {£\/ kH)"^ is in 27rN and large, then, 

c{eV^) = fKH{E2 + o(l)) . 
Inserting this into ()6.12p and using that Me ~ we get. 



£o{hi;,A) > {1-6)\D\[k- H]l{{l + 2a)E2-Ca~^{K- H) 



2n2\ 



IV'pdx. (6.13) 



Also, using Holder's inequality and Lemma 16.21 we get ||V'lli2(-f7) ^ CQ- That way, we infer from 
([UT^ the bound, 

Eo{hij,A) > {1 - 6)\D\[k - H]l {{1 + 2a)E2 - Ca-^^C'^f) 

-c((5 + rtC^/V + ^)c. (6.14) 

Let 

ei = max(|ln(C"/^K)|-i,C^/8). 



3D GINZBURG-LANDAU FUNCTIONAL 



33 



It results from the definition of the parameter that ei ^ 1 as k — )• oo. We select the parameters 
S, a and i as follows, 

a = ei, i= 7=, 6 = eiC- 

It is easy to check that, 

K^C^^^ < 1 , ^ + 5«^C + S'^^C^^^K^ < , and < £ , 

as K oo. 

We select the parameter a as follows, 



-3,^2a2^2 



so that fj <C 1 and a k % £ ^ 1 as k — t- oo. Thanks to this choice of parameters, and noticing 



that K^C'^ = max(ft;, [k — H]'^), we infer from ()6.14p . 



£o{hi;,A) > \D\[k- H]lE2 + o(^max{K, [k - H]l)^ 
This finishes the proof of Theorem 16.11 



□ 



6.2. Energy estimate of a trial configuration. 

The aim of this section is to prove Theorem 16.51 below, which estimates the energy of a test 
configuration. The construction of the test configuration requires some notation. Let x be a 
cut-off function such that, 

XGC^m, suppxC [-2,2], 0<x<l inM, X = l m[-l,l]. 
Let T] G (0, 1) and define a function /i^ G C^{il.) as follows, 

/ dist(x, (9il) 



V X G , h„(x) = 1 — X { - 

V V 

Notice that vanishes in a thin neighborhood of the boundary, 

Q(r]) = {x edn : dist(2;,a0) < r?} , 

and = 1 in i} \ n{2ri). 

Consider a parameter £ G (0, 1). Let b = H/k and R = Recall the ground state energy 

mp(b,R) defined in (|2.21|) . together with the space Eji in (|2.2Up . Let u^^ji G Eji be a minimizer 
of the energy in (|2.ip . i.e. 



mp{b,R) = ^6|(V - iAo)ufe^ijp - |ub,Rp + ^lufe^fll^^ dx , 

where Kr = {-R/2,R/2) x {-R/2,R/2) and Ao(xi,X2) = (-X2/2, xi/2). By definition of the 
space Er, the function Ub,R G i7j^^(]R^). 

We define the test function V'Jj'^ G H'^{Q;C) as follows, 

yxen, i^^%{x) = h^{x)ub,R{^x±). (6.15) 

Here we used the notation that if 2; = {xi,X2,X3) G M^, then x± = {xi,X2)- 

Theorem 6.5. Suppose that the magnetic field H is a function of k such that, 

H H 
1 < lim inf — < lim sup — < 00 . 

Suppose furthermore that I and ij are functions of k such that, as k ^ 00, 
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Let R = £^/kH, be the function introduced in ()6.15p and F{x) = (— X2/2, xi/2, 0). Then, 



as K ^ 00, 



£{'^^%,F) < E2\n\ [k - H]l + o(^max(K, [k - H]l)^ . 
Here E2 < is the universal constant introduced in (j2.42p . and E is the functional in p.ip . 

Proof. Recall the magnetic potential Ao(xi,X2) = {—X2/2, xi/2) and that the function Ub^R 
satisfies the equation , 

-&(V - iAo)ub^R = (1 - \ub,R\'^)ub,R in . 

1/2 

By assumption on the magnetic field, we can make [1 — arbitrarily close to for all k > K^hy 
choosing kq > sufficiently large. It results from Proposition 12. 6l that u^^r satisfies the universal 

bound \ub,R\ < Cmax[l - b]^/'^. 
An integration by parts gives us, 

/ {b\{V-iF)fub,Rf-\fub,R\^ + f^\ub,R\^)dx=[ \\Vf\ub,Rfdx, (6.16) 



for any compactly supported smooth real- valued function / G C|(M'^) satisfying < / < 1 in 
M.^. Consequently, we get 

/ (b\(y-iF)fub,R\^-\fub,R\'' + l\fub,R\A dx<l [ f\f^ - 2)\ub,R\Ux 

+ CLx[l-b]+ / iV/pdx. (6.17) 

Let D = supp/ and {QR,j)j£z'^ a uniform lattice of M.^ whose unit cell is the cube 

Qr = Krx {-R/2, R/2) , Kr = {-R/2, R/2) x (-i?/2, R/2) . 

Let J = {j : Qrj n I? / 0} and N = card J. Recah that \ub^R\ is periodic in Then, 

writing 

/ f\2- f^)\ub,R\Ux= f \Ub,R\Ux+ [ {2f-f^-l)\ub,R\^dx 

= Nr[ \ub,R\Ux^+ [ {2f^ - f^ - l)\ub,R\Ux . 

1/2 

We use the bounds |u6,_r| < Cinax[l — &]+ and < / < 1, and we notice that an integration by 
parts gives mp{b,R) = —'^fj^^ \'Ub,R\'^dx. Consequently, we obtain that, 

/ f\2-f^)\ub,R\Ux>-2NRmp{b,R) + 3Ci^,[l-b]l[ {f^-l)dx. 

Inserting this into ()6.17p we get, 

j^^ (b\{V - iF)fUb,R\^ - + ]^\fUb,R\^^ dx 

<NRm,{b,R) + CLAl-bh [ [^-^]+ (l-/2)) dx. (6.18) 
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We choose f{x) = /i^(x/V kH). Performing a change of variable in (j6.18p . we get (recah that 
b = H/ K by definition), 

H\/kH f f rrT7\„/,blk |2 ,.2|„/,blk|2 , '^^|,,blk|4 



'ri,R\ I dx 



< NRm,{b, R) + MCLAl - bh + !I!lE^El^LJ\±^ (g.lQ) 

where M is a universal constant. With this choice of /, it is easy to check that the number 
satisfies, 

X ( I — > ir^l , as i? — > oo . 



We insert this into (|6.19p . Theorems 12.111 and 12.121 together tell us that 

mp{b,R) < R^[l-b]l{E2 + o{l)) as R ^ oo , 
where E2 is a universal constant. Therefore, we deduce from (|6.19p that, 

£:«iF) < .2[1 - b]lE2il + o(l)) + M'CLx[l -b]J^ + ^^maxK^il - bl 



r] 2 



where M' is a new universal constant. 

Since k"^ ^ ^ 1 as k — )• 00, it is easy to check that 

t^«A^2maxfi,[l-6]2 

Since b = H/k, this finishes the proof of the theorem. □ 
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